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Abstract 



The natural relativistic generalisation of Landau's two constituent superfluid 
00 ■ theory can be formulated in terms of a Lagrangian C that is given as a function 

g; -^^^^^-^ 

■ phase scalar. It is shown that in the "cool" regime, for which the entropy is 
attributable just to phonons (not rotons), the Lagrangian function C{s, V92} is 
given by an expression of the form C = P—3ip where P represents the pressure 

■ as a function just of V p 92 in the (isotropic) cold limit. The entropy current 



dependent contribution ip represents the generalised pressure of the (non- 



■ isotropic) phonon gas, which is obtained as the negative of the corresponding 

grand potential energy per unit volume, whose explicit form has a simple 
algebraic dependence on the sound or "phonon" speed c p that is determined 
by the cold pressure function P. 
PACS numbers: 04.40, 67.90 

I. INTRODUCTION 

The purpose of the present work is to derive the natural cool limit form of the equation of 
state that is needed to complete the formulation of the natural relativistic generalisation of 
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Landau's two constituent superfluid theory J]J . The qualification "cool" is to be understood 
here as referring to the low temperature limit in which account is taken only of lowest order 
deviations from a nearby cold configuration in which thermal effects are absent altogether. 

Whereas the cold (strictly zero temperature) superfluid case is described by irrotational 
configurations of a simple perfect fluid model (i.e. one that is both isotropic and barytropic), 
on the other hand the allowance for thermal effects, even in the cool limit, necessitates 
the use of a two constituent fluid model, of the kind that was pioneered by workers such 
as London, and perfected by Landau. Accurate treatment of non-stationary applications 
would require allowance for viscosity of the "normal" constituent. However it is sufficient 
for many purposes - and entails no loss of accuracy at all in the case of equilibrium states 
such as simple cylindrical vortex configurations - to use a strictly conservative treatment as 
in Landau's original model Jl|. 

In addition to the neglect of dissipation, the main simplification on which the analysis 
below is based is the neglect of the non-linear excitations known as "rotons" which are of 
dominant importance in the "warm" regime nearer to the phase transition (beyond which lies 
the regime of "hot" states in which superfluidity is absent altogether). The presence of such 
effects makes it very difficult to derive the equation of state over the full range of the three 
relevant variables, which in addition to the temperature G say, could be taken to consist of 
the relevant conserved particle number density n say, and the relative velocity v of the two 
constituents. Although it is out of the question for neutron star matter, an experimental 
investigation of the equation of state should be feasible in practice for the laboratory case 
of Helium-4. However even in this experimentally accessible case, the exploration of the 
relevant phase space (as parametrised by 0, n, and v) does not yet appear to have been 
sufficiently systematic and thorough (see e.g. 0). 

In contrast with the dubious quality of present day knowledge of "warm" superfluid 
dynamics, the state of knowledge of the "cool" limit is very satisfactory. In this limit, 
complications such as "rotons" can be ignored, the only important thermal effects being 
entirely attributable to simple "phonon" excitations, which can be adequately described by 
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an essentially linear treatment whose original development is again largely attributable to 
the (by now experimentally well substantiated) work of Landau ||. The detailed devel- 
opment of the necessary statistical mechanics has been summarised in a convenient form 
by Khalatnikov |4[]. It will be shown below that although it was originally carried out in 
a non-relativistic framework, the nature of this statistical analysis is such that it can be 
translated into a relativistic form without any change of form provided that sufficient care 
is used in defining the appropriate variables for the relativistic version. In consequence, as 
shown in the appendix, Landau's expressions for the first and second sound speed can also 
be retained without change in the relativistic regime. 

II. RELATIVISTIC GENERALISATION OF 2-CONSTITUENT SUPERFLUID 

DYNAMICS 

It has recently been shown [4] that the natural relativistic generalisation of the standard 
Landau theory of two constituent superfluid dynamics in its simplest strictly conservative 
version [1] can be formulated in a particularly convenient manner by taking as the starting 
point a Lagrangian scalar £ that is given as an appropriate function (which was denoted by 
X in the original presentation ||) of the entropy current vector s M and of the gradient 



of the superfluid phase scalar (p. Two of the required equations of motion are obtained in 
the form of the usual conservation laws as given by 



(2.1) 



(2.2) 



for the particle current n p constructed below, and 



V p s p = . 



(2.3) 



The system is completed by the equation 
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s p V [p e a] = (2.4) 

(using square brackets to denote index antisymmetrisation) governing the evolution of a 
certain thermal momentum covector Q p . This last equation can alternatively be expressed 
(in a form that is useful as a starting point for the derivation of corresponding circulation 
and helicity conservation laws ||) as the vanishing of the Lie derivative of the thermal 
momentum vector with respect to the associated temperature vector, j3 p , i.e. 

(3 a V a Q p + e a V p (3 a = , (2.5) 

where 

f3 p = -( s ff e ff rV . (2.6) 

The thermal 4-momentum covector p and the particle number current vector n p are spec- 
ified in this formulation as the algebraic functions of the entropy current vector s p and of 
the superfluid 4-momentum covector fi p that are obtained by partial differentiation of the 
Lagrangian according to the infinitesimal variation formula 

alC = & p ds p - n p d/j p . (2.7) 



III. ALTERNATIVE ACTION FUNCTIONS AND THE STRESS MOMENTUM 

ENERGY TENSOR 

The formulation that has just been summarised has the technical advantage of being 
particularly economical in so much as it involves only 5 independent component variables, 
namely the phase scalar (p and the 4 independent components of the entropy current vector 
s p . This feature of economy has recently been exploited for the purpose of setting up a 
correspondingly economical Hamiltonian formulation of the theory J7|, and it has also been 
exploited as a guide to the formulation of an analogously economical theory for describing 
thermal effects in superconducting cosmic strings f|. 
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Although mathematically equivalent as far as its physical ( "on shell" ) solutions are con- 
cerned, this more recent 5-component formulation differs from the earlier convective varia- 
tional formulation P,|10f (which has the alternative advantage of being less specialised) that 
treated the current vectors n p and s p on the same footing as independent variables, thus 
involving a total of 8 independent space time components. The presence of the three extra 
components made it necessary to include, as an extra dynamical equation, the irrotationality 
condition 

V [p /V] = (3.1) 
that is the Poincare integrability condition for the existence of a potential tp satisfying 



(2.1). In the convective variational formulation jj, p is on the same footing as B p , both being 



obtained by partial differentiation according to the formula 

dA = Q p ds p + fi p dn p , (3.2) 

where the convective variational master function A is obtainable from the Lagrangian £ of 
the newer more specifically adapted approach by a Legendre type transformation given by 
the relation 

A = C + r//i P - (3.3) 

The approach based on the Lagrangian £ that will be used here is in fact a compromise, 
intermediate between the convective variational approach based on the master function 
A, and another independently developped (Clebsch type) variational approach |TT| , |12"11 in 
which both currents are treated as dependent variables, their specification being given by a 
generalised pressure function \l/ in terms of as defined by ( |2.1| ) and of O p , which in this 
version has independent status, by the partial differentiation formula 

dty = -s p dQ p - n p dfi p (3.4) 

where the generalised pressure function \1/ itself is given by 



V = C- e p s p . (3.5) 

Although the generalised pressure function \l> will not act in a fundamental role in the 
formulation used here, it nevertheless plays an important part, notably in the expression 
for the stress momentum energy tensor. An important part is also played by a set of three 
scalar functions of state, C, B, and A (respectively describable as the caloric coefficient, 
the bulk coefficient, and the anomaly coefficient), that are defined as matrix elements in the 
relation expressing the momenta in terms of the velocities in the form 

Q p = Cs p + An p , // = As p + Bn p . (3.6) 

These coefficients can be used to convert the canonical expression 

T p a = n p fi a + s p Q a + Vg p a (3.7) 

for the stress momentum energy density tensor of the 2-constituent superfluid into the 
manifestly symmetric though not so elegant form 

T P<J = Bn p n a + A(n p s a + s p rf) + Cs p s p + ^g pa , (3.8) 

where g pa is the contravariant inverse of the - flat or curved - space time metric tensor g pa 
that is to be used for index raising or lowering. For the purpose of the present work it will 
be more convenient to work with an alternative recombination expressible in terms of the 
dilation coefficient $ 2 and the determinant coefficient K, by 

T pa = ^(//y 7 + £s p s a ) + ^g pa , (3.9) 

where 

& = B- 1 , K = CB-A 2 . (3.10) 
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IV. THE COLD LIMIT AND THE SPEED OF SOUND 



Before considering the treatment of the lowest order thermal corrections, it is to be 
recalled that in the cold (i.e. strictly zero temperature) limit, the superfluid behaves as a 
perfect fluid of simple barytropic type. This means that in terms of the unit flow vector u p 
defined by 

n p = nu p , u p u p = -c 2 , (4.1) 

where c is the speed of light, the stress tensor reduces to the isotropic form 

T pa = (p + c- 2 P)u p u a + Pg pa , (4.2) 

in which the mass-density p and the (isotropic) pressure P are directly related by a single 
variable equation of state. The equation of state can be specified either by giving the rest 
frame energy density p as a function of the corresponding particle number density n in the 
form 

p = p{n} , (4.3) 

or equivalently, in dual form, by giving the pressure P as a function of the corresponding 
effective mass (i.e. relativistic chemical potential) variable p in the form 

P = P{p} , (4.4) 

(using small curly brackets to indicate functional dependence, as distinct from the multipli- 
cation that would be indicated by ordinary brackets). The connection between these two 
formulations is given by the familiar differential specifications 

p — — n- — — (4 5) 

dn ' c 2 dp ' 

and the dually symmetric relation 

p + c~ 2 P = np . (4.6) 



For the simple ( "barytropic" ) perfect fluid model that has just been formulated, the 
convective variational master function A referred to in the introduction is given directly by 
the first version of the equation of state as 

A = - c 2 p{n} (4.7) 

while the corresponding dual Lagrangian C, which in this zero temperature limit has just 
the same form as the Clebsch type potential function will be given analogously by the 
dual version of the equation of state as 

C = V = P{/j} . (4.8) 

To complete the translation into the notation of the preceeding section the particle momen- 
tum covector will be given simply by 

Hp = (M p , (4.9) 

and of course the entropy vector s p and the thermal momentum covector P will both vanish, 
while the temperature vector f3 p is singular in this zero temperature limit. 

One of the most important quantities in this simple perfect fluid model is the charac- 
teristic sound or "phonon" speed, c p say, i.e. the short wavelength limit of the propagation 
velocity relative to the preferred rest frame of small (necessarily longitudinal) perturbations. 
This characteristic speed is immediately derivable from either of the versions ( |4.3| ) or ( |4.4| ) 
of the equation of state via the familiar differential formula 

c 2 = d 4- = c 2 ^ . (4.10) 
p dp nan 

This can be used to construct the phonic (or sonic) metric tensor 

p * = *r + (^-^)t*v r , (4.11) 

whose null eigen co vectors p p , as defined by 

Q pa p P Pa = , (4.12) 



are tangential to the characteristic hypersurfaces of sound propagation in the medium, in 
the same way as the null covectors of the ordinary spacetime metric g pa are tangential to 
are the characteristic hypersurfaces of light propagation in vacuum. 

V. THE CONCEPT OF "NORMAL" AND "SUPERFLUID" DENSITY 

CONTRIBUTIONS 

As soon as we want to allow for deviations from the cold limit case that has just been 
summarised, it is necessary to increase the number of independent variables in the funda- 
mental equation of state from 1 to 3. It might naively have been hoped that 2 would suffice 
for the cool limit with which we are concerned here but, as will be made apparent at once, 
this ceases to be possible unless one is willing to restrict attention to states for which there is 
no relative motion between the two constituents. In the convective variational formulation 
based on A as the fundamental state function, the 3 independent variables correspond to 
the 3 scalar invariants that can be constructed from the pair of independent current vectors, 
namely n p whose direction determines the particle or "Eckart" rest frame, and s p which 
similarly determines the thermal or "normal" rest frame, so that an obviously convenient 
way of choosing the 3 independent variables in a fundamental state function of the form 



will be to take them to consist of the particle rest frame number density n as given by 



A = A{n, x, s} 



(5.1) 



c 2 n 2 = — n p n p , 



(5.2) 



the cross product variable x as given by 



c 2 x 2 = —n p s p , 



(5.3) 



and the thermal rest frame entropy density s as given by 



In the dual formulation based on the lagrangian £ that will be used here, the relevant 
independent variables are the 3 scalar invariants that can be constructed from the entropy 
current vector s p and the momentum covector fi p whose orientation determines what is 
known in this context as the "superffuid" rest frame. It is obviously convenient to take one 
of the required scalars to be again the "normal" entropy density variable s as defined by 
( |5.4|) , while taking the other two to consist of a new cross product variable y given by 



cV = -fi p s p , (5.5) 

together with the effective mass variable fi given by 

c 2 /i 2 = — /i p /i p , (5.6) 

so that the required generalisation of ( [4.8D , i.e. the relevant analogue of (|5.1|) , will be given 
by an expression of the form 

£ = £{v,y,s}- (5.7) 

Starting from an expression of this form it can be seen from ( |3.6p and ( |3.10| ) that the 
secondary variables n p and P will be given in terms of the primary variables fi p and s p of 
this formulation by 



n 



p 



$ 2 (// - As p ) , 9 P = $ 2 (/Cs p + Alio) • (5.8) 



where by ( |2.7|) the relevant dilation, determinant and anomaly coefficients $ 2 , /C, and A are 
given by the partial differentiation formulae 

c 2 $ 2 = -^, c 2 $ 2 /C = -^, c 2 $M = -l^. (5.9) 
/i on s os 2y oy 



It was a generic Lagrangian function of this form ( |5.7| ) that was the starting point for the 
construction, in implicit form, of a Hamiltonian reformulation of the theory 0. However in 
order to be able to make such a construction explicit it is necessary to sacrifice generality and 
deal separately with particular kinds of state function. From a purely mathematical point of 
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view the simplest kind that can be envisaged is the "regular" category, as characterised JTG 
by the condition that the anomaly parameter A vanishes, which evidently, by the original 
definition ( |3.6| ), means that the momenta are aligned with the corresponding current vectors. 
It can be seen from ( |5.9| ) that in the formulation based on the Lagrangian £ that we are using 
here this regularity condition is expressible as the condition that the state function should 
depend only on p and s but not on the third cross product variable y. However although 
it can be hoped that such a mathematically attractive alignment ansatz may provide a 
useful approximation in contexts involving hot conducting fluids, it is certainly not at all an 
appropriate simplification in the context of superfluidity. 

In order to see how to obtain a more suitable form of equation of state it will be useful 
to establish an appropriate translation relating the terminology of the traditional formalism 
that was developed in a non-relativistic framework to that of the relativistic formalism used 
here. In particular it will be convenient to introduce appropriate relativistic generalisations 
of the concepts of the "normal" and "superfluid" mass densities p N and p s . One way of 
defining such a generalisation is to base it on a corresponding vectorial decomposition |5j 
of the particle current n u , which in the non-relativistic limit is interpretable as being 
proportional to the Newtonian mass current p u = mn u for some constant coefficient m 
representing the mass per particle. However such a definition has the disadvantage that in a 
relativistic context there will be a certain margin of ambiguity in the choice of the appropriate 
"rest mass" parameter m. It will be more convenient for most purposes, including that of 
the present work, to adopt a natural alternative possibility that gives the same result in the 
non-relativistic limit but that is free of any such ambiguity of normalisation. 

The definition that will be adopted here - and that seems clearly most appropriate in the 
relativistic context - is based on the decomposition of the stress momentum energy density 
tensor in the form 

^ = ^7 + ^/ + $^ , (5.10) 

/i N s s 

where s s is the value of the entropy density in the "superfluid" rest frame defined by p p , 
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and /i N is the analogously denned superfluid chemical potential with respect to the "normal" 
rest frame defined by s p , i.e. 

/, v\-i/2 y 2 / v\-i/2 y 2 

s s = s \ 1 --i) ^ = H 1_ ^J =7' ( 5 - n ) 

in which the velocity t> is the relative translation speed between the "normal" and "super- 
fluid" frames as given by 



The definitions implicit in (|5.10 ) are such as to ensure that the momentum density (T 



say) in the direction of the relative motion will have magnitude p N v in the "superfluid" 
frame, while it will have magnitude p s v in the in the "normal" frame. (One could define a 
corresponding generalisation of the Newtonian concept of the mass current by setting p u = 
Ps/^n -1 / 1 ^ + Pn s s~ 1sU > ft wou ld be of limited utility since it would not retain the exact 
conservation property of its Newtonian limit.) It can be seen by comparison with (|3.9[) that 



the values of the quantities so defined can be evaluated in terms of the partial derivatives 
appearing in (|5.9| ) using the formulae 

p s = $ 2 /i N 2 , P N = $ 2 /Cs s 2 . (5.13) 



VI. STATISTICAL MECHANICS OF THE PHONON GAS 

As in the standard Landau analysis 0, it will be supposed that, in the cool limit with 
which we are concerned, the thermal excitations can be described as a bosonic gas of non 
interacting phonons. 

To see how this works out, and to fix the notation, we start by recalling that the esssential 
step in such an analysis || is the evaluation of the expected value of the occupation numbers 
u r , say where r is an index running over the relevant quantum states. The standard trick for 
doing this is work in terms of groups, with collective index R{r} say, consisting of a large 
number, G R say, of nearby states characterised by approximately equal quantum numbers, 
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the point of this being to ensure that although the mean occupation number v R may be 
small, the total group occupation number, 

K = G R u R = £ u r (6.1) 

r:R{r}=R 

will be large compared with unity. This makes it possible to replace the exact Bose statistical 
formula 

for the total entropy S of the whole system (in units such that Boltzmann's constant is set 
to unity) by the corresponding Stirling formula 

S = E G r (fa + 1) ln K + 1} - ln **) • ( 6 - 3 ) 

R 

Introducing the dimensionless state functions cr{r} = a{R{r}} defined by setting 

1 

'« _ e *{R} - 1 . u 



v r , = ■ ,„| , . o {/.'} = ln {^-^} , (6.4) 



R 



the formula (6.3) can be usefully rewritten as 

S = lnZ + ]TA/>{i?} , (6.5) 
/? 

where the zustandssumme, Z, is defined as a sum over all admissible combinations {u r } of 
individual occupation numbers by 

z = ex p{XXo-{r}} , (6.6) 

K} 

from which the logarithm required in fl6.5| ) is obtained simply as 

lnZ = £G fl ln{^ + l}. (6.7) 
The corresponding variation formulae are thus obtained even more simply as 

5\nZ = ~Y,vM{r} , = ^>{r} 5z/ r . (6.8) 
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Having carried out these routine preliminary steps, one can immediately obtain the 
thermal equilibrium state specified by maximising the entropy subject to possible constraints 
on the total occupation number N and the total energy momentum covector V p as defined 
in terms of the corresponding microscopic energy momentum covectors p p by 

JV=£^* = £",., V P = J2vrP P {r}- (6-9) 

R r r 

In order to satisfy the ensuing requirement 

5S + a8H + P P 8V P = , (6.10) 

where the scalar a and the 4 vectorial components (3 P are Lagrange multipliers, it evidently 
follows from ( |6.8|) that the solution for the equilibrium distribution must be given simply by 

a = -a - p p p p , (6.11) 

which by ( |6.5| ) implies the thermodynamic relation 

S-\nZ = -aU-l3 p V p . (6.12) 

In the case of phonons, as in that of the more familiar example of Planckian photons, 
there is no conservation law imposing any constraint on the total occupation number so the 
multiplier a will vanish, while the vector (3 P will give the temperature and relative flow 
velocity components v i (i— 1,2,3) with respect to the chosen reference frame according to 
the specifications 

a = 0, P° = ^ P^y*- (6 ' 13) 

For the actual evaluation of the distribution thus obtained it is of course convenient to 
work in the continuum limit for which the summation goes over to a phase space integration 
in the simple form 

E^ = E^V/^, V = Jdx^dx*. (6.14) 



R 



(2vrft) 
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This is the first point at which the phonon case with which we are concerned here deviates 
from the more familiar photon gas, in which an extra factor of 2 would be required to allow 
for the existence of 2 (transverse) polarisation modes, whereas in the phonon case there is 
only a single (longitudinal) polarisation mode. 

The only integration for which ( |6.14| ) is actually needed in practice is that for the loga- 
rithm of the zustandssumme as given by (|6.7f ), since once this is known the other thermody- 
namic quantities involved will be obtainable from it by straightforward partial differentiation 
procedures. Defining the grand potential energy i?, the free energy F, and the ordinary en- 
ergy U of the system by formulae of the standard form [|3| 

Q = -Q\nZ, F = U-QS, U=-V , (6.15) 

(which differs however from the terminology of Khalatnikov |4[] who uses the symbol F in 
place of i?) the relation (6.12) can be rewritten as 

F-Q = v i V i , (6.16) 

while the variation laws ( |6.8p (which are of course to be understood as determined with 
respect to fixed boundary conditions, and hence in particular for a constant value of the 
volume V) give rise to corresponding variation formulae of the form 

5H = S8e -Pi 8v i , 5F = -SSQ + v^Vi , 5U = G8S + v i 8P i . (6.17) 

Although it is simpler in so far as polarisation is concerned, on the other hand the phonon 
case is more complicated than that of photons in that the presence of the background medium 
implies a breakdown of Lorentz invariance which is expressed by the need to use not the 
ordinary metric but the phonic metric ( [4. 1 1| ) in the relevant nullity restiction ( |4.12| ). In 
order to be able to satisfy this condition by expressing the phonon energy e by the simple 
formula 

e = -p = c p p , p 2 = pi + p 2 2 + pi (6.18) 
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where c p is the zero temperature sound speed as given by ( |4.10|) , we must now restrict 
the choice of reference frame to be one in which the underlying superfluid is at rest. This 
contrasts with the case of a photon gas, which is subject to an analogous formula, with 
c replacing c p , in an arbitrarily boosted frame. Further restricting the frame by choos- 
ing space axes aligned with the relative flow direction, one obtains the logarithm of the 
zustandssumme, in the explicit form 

, „ f 2txv 2 sin 6 d9 dp , . 

ln ^- v / M1 " e } ' (619) 

with 

a = = £^ , (6 . 2 o) 

where 

, v 3 = v , p 3 = pcos9 . (6.21) 



VII. THERMODYNAMICS OF THE PHONON GAS 

To relate the thermodynamic quantities obtained in the preceeding section to the su- 
perfluid continuum variables described in the earlier sections, we start by introducing the 
generalised pressure function ip determined by the grand potential fi of the phonon gas, 
together with the entropy density s s in the superfluid frame - which is of course to be 
identified with the component given by ( |5.11| ) - according to the specifications 

n = -Vtfj , S = Vs s (7.1) 

while the corresponding preferred frame components £ and II of the phonon energy and 
momentum density are specified by setting 

V = -VS , V :i =Vn , V 2 = V, = . (7.2) 



With these definitions the global thermodynamic equilibrium condition ( |6.12| ) can be rewrit- 
ten in terms of the corresponding local field variables as 
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£ + ip = es s + nv, (7.3) 

while by ( |6.17| ) the fundamental entropy maximisation condition ( |S.1(J| ) gives the relevant 
local version of the first law of thermodynamics in the form 

58 = e5s s +v5n , (7.4) 

whose conjugate, by ( f7T3|) , is evidently 

Sip = s a 5Q + n5v . (7.5) 

Using this last formula one can obtain the values of s s and II as functions of and v by 
partial differentiation of the generalised pressure function if), which is obtained from (|6.19|) 
in the form 

(7 - 6 » 

where h is a constant that is given, within a numerical factor that is extremely close to unity, 

by the usual Dirac Planck constant h, its exact expression being 

~ / 1215 xi/3 t 9 /5ttx 1/3 ^ ^ . , 

h=( -) h= — (— ) ft- 0.99ft . 7.7 

One thus obtains the expressions 



s 



and 



which are formally identical with the analogous expressions as originally derived in a non- 
relativistic framework 0]. 

It is to be remarked that the temperature vector (|6.13 ) can be used to rewrite ( |7.6|) in 
the neater form 

^ = (^) 3 c p (2^^)- 2 , (7.10) 
where Q~p a is the covariant inverse of the sonic metric given by ( |4.11| ). 
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VIII. CONCLUSION: THE EQUATION OF STATE FOR COOL SUPERFLUID 

DYNAMICS 



We can use the expressions obtained in the previous section to evaluate the "normal" 
mass density contribution p N as introduced in (|5.10|) in a manner that is formally consistent 



with the traditional definition || by identifying the momentum density II with the corre- 
sponding mixed component of the stress momentum energy density tensor in the preferred 
superfluid frame defined by fi p according to the specification 

T° 3 =II = p N v. (8.1) 

This leads to a result that is expressible in the form 

We are now in a position to obtain the equation of state for the required cool superfluid 
Lagrangian £ by integrating the partial differential equation 

ldC c 2 



s ds s c 



2Pn (8-3) 



obtained from ( |5.9| ) and ( |5.13| ). This leads to the principal result of this work which is the 
formula 

C = P{fi}-3^{/i,y,s} , (8.4) 

where ip is the generalised pressure function of the phonon gas, which can be seen from 
( pM0| ) to be given by 

3V> = fic-^g-^s*^ . (8.5) 

Since the operation d/ds is defined in terms of the scalar variables /i, y, and s, the verification 
that ( |3.4| ) satisfies ( |3.3| ) requires the conversion of ( |S.5| ) into terms of these three scalars. 
This can be done using the explicit expression 
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1 c 2 

= 9p° + ^(l - ~^)UpUa , (8.6) 



for the covariant version of the phonic metric ( |4.11| ), which gives 



e> ff = (c 2 -c:)(^-cV. (8.7) 

The solution of ( |3.3|) is of course not unique (since one can always add in any function of 
fx and y) but ( |8.4j ) is the only solution compatible with the boundary requirements that it 
should go over correctly to the cold limit value ( |4.8j ) when both s and y vanish, and that the 
energy in the static limit v — should differ from that of the cold solution with the same 
particle number density n by an amount given by the phonon gas contribution £ , which in 
the static limit is the same as 3i/j. 

We conclude by remarking that the description of ip as the generalized pressure function 
of the phonon gas is justified by the fact that the basic cool equation of state formula ( |S.4|) 
translates into terms of the total generalised pressure function as defined by ( [3751) , simply 
as 

* = P + V, (8.8) 
where the appropriate expression for the thermal correction term ip is 

^ = T(|) 3 (^ ^) 2 ' (8-9) 

in which it is to be recalled that (as described in Section 4) c p and Q pa are given as functions 
(only) of the superfluid momentum covector fi p in a manner determined by the original 
equation of state for the pressure P in the zero temperature limit. The corresponding 
expression for the relationship between the entropy current vector s p , the temperature vector 
j3 p , and the thermal momentum covector P is 

s p = A^(3 p = 2(i) 3/2 (c p ^) 1/ V CT (7 . (8.10) 

The first of these relations can be used to verify that the temperature vector (3 P introduced 
as a Lagrange multiplier in Section 6 agrees with the earlier definition (|2.6j ). 
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APPENDIX: EVALUATION OF THE FIRST AND SECOND SOUND SPEEDS 

To obtain the characteristic speeds of propagation of small perturbations, it suffices to 
follow the lines already developed in the context of more general multiconstituent fluid the- 



ory [10] using a technique originally due to Hadamard. The speeds in question are those 
characterising the set of possible directions for the normal covector, X p say, of a charac- 
teristic hypersurface of the system. The Hadamard analysis postulates that the relevant 
independent dynamical variables, in this case s p and fi p , should be continuous across the 
characteristic hypersurface, but that their derivatives have infinitesimal discontinuities that 
will be expressible in the form 

pv> CT ] = \ p s a , [V p /v] = \ P K , (Al) 

in terms of a corresponding set of infinitesimal discontinuity amplitudes s p and ft p . For the 
dependent variables n p and P , one then obtains the analogous relations 

[V p n CT ] = \ff , [V P 6 CT ] = A P CT . (A2) 

For a generic Lagrangian of the form (|5.7|) the independent dynamical variables fi p and 
s p will determine the dependent dynamical variables n p and G p by the relation (|5.8|) , which 
can be rewritten succinctly as 

n p = Bfi p - As p , Q p = Cs p + A(i p , , (A3) 



with 



B = 2 9C , C = -2-^- , A = 9C (A4) 
c 2 <9/i 2 ' c 2 ds 2 ' c 2 dy 2 
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It follows that the independent discontinuity amplitudes jj p and s p will determine the cor- 
responding dependent discontinuity amplitudes h p and P by an analogous linear relation 
of the form 

rip = Bpufi — A ap s , Op = C pa s + Apufi , (A5) 

with 

dB dA dA 

Bpa — Bg pa - 2^ fipfig. + 4 c2 ^ 2 /i(pS CT ) + c2 ^ 2 s p s a , 



dC dA dA 

^ pa - ^ 9 pa Z C 2Q S 2 S P S(T 4 c 2Q g 2 S (P^) c 2Qy2^P^ a ' ^ At V 



_ ) _8A_ _ dA _ dC _ dC 

A pa - A 9pa ^/V* ^g^ S P S - 2 JQ^ Sp ^ 



The required characteristic equations are to be obtained by substituting the formulae 
J]) and ( |A2| ) in the discontinuities of the relevant dynamical equations, which are the 
integrability condition ( J3.1Q for together with (|2.2|) , ( |2.^| ) and (P^). The last of these 



gives 

s p \ [p Q a] = , (A7) 
the flux conservation equations fl2.2|) and (|2.3| ) give 

A p n p = , (A8) 

and 

\ p s p = , (A9) 
while finally the irrotationality condition ( |3.1|) gives 

\[ P K\ = , (A10) 

which is interpretable as meaning that fi p must be proportional to \ p . The weaker condition 
gives rise to two qualitatively different possibilities, the one of interest in relation to 
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the first and second sound modes under consideration here being that P should also be 
proportional to X p . This condition is expressible by replacing (|A7|) by the formal analogue 
of ( [MO) ), i.e. 

A[A] = • (All) 

The alternative possibility is the trivial case of a thermal shearing mode as characterised 
by s p \ p = and s p Q p = 0. In the more general non-superfluid case []10| the ordinary 
particle constituent could also have a shear type mode, but this second kind of trivial mode 
is excluded in the superfluid case by the irrotationality condition ( |A10| ). Leaving aside this 
trivial case, we now restrict our attention to the sound type modes characterised by the 
restriction of flA7|) to flATTl ). 

To deal with the ensuing system of characteristic equations (|Ai|), (|A9|) , (|A10|) and ( |A11|) , 



we now put it in a more explicit form by the use of local "superfluid frame" Minkowski 
coordinates {x p } = {x°, x 1 , x 2 , x 3 } aligned with the relative flow direction at a particular 
point under consideration, so that the independent dynamical variables will be given there 
by 

{//} = /.{I, 0, 0, 0} , {s p } = s g {l, v, 0, 0} . (A12) 

There will be no loss of generality in taking the characteristic covector to have the form 

{X p } = {—u, cos#, sin#, 0} , (A13) 

where 9 is the angle between the relative flow direction and the direction of propagation 
of the perturbation, and u is the characteristic velocity whose evaluation is the ultimate 
objective of the exercise. The ensuing set of characteristic equations can be organised in 
three subsets. Firstly there is a "superfluid" subset obtained from (|A8|) and (|A10|) in the 
form 

uh° — cos On 1 — sin 9 n 2 = , (A14) 
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cos 8 ji + u $i x = , sin 9 fi 1 — cos 9 fi 2 = . (A15) 
Next there is a "thermal" (or "normal") subset obtained from ( |A9| ) and ( |A11| ) in the form 

us — cos 9 s 1 — sin 9 s 2 = , (A16) 



cos#9 = , sin^6 1 -cos#6> 2 = . (A17) 



Finally there is an orthogonal subset containing the remaining information from ( |A10|) and 
( |A11|) in the form 

/2 3 = , 6 3 = . (A18) 

This last set ( [A18Q will always decouple and can be solved separately so as to provide 
the conclusion that orthogonally to the plane defined by the relative flow and propagation 
directions, not just the momentum discontinuty amplitudes but also those of the currents 
must vanish: n 3 = s 3 = . 

In the limit when the temperature G and the entropy magnitude, s = 0{<d 3 } go to zero, 
one is left only with the first subset of characteristic equations (|A14|) , ( |A15|) corresponding 
to ordinary sound modes, but when entropy is present it will evidently be necessary to take 
account also of the second "normal" subset (|A16| ), ( |A17|) . In general this second set will be 



coupled with the first by the cross terms in flAq), so that one will obtain a rather intractable 
quartic characteristic equation for the propagation speed u. However it transpires that in the 
"cool" limit with which the present work is concerned, i.e. to lowest order in the temperature 
O, the first and second subsets will decouple, giving a pair of quadratic equations that can 
easily be solved to give the respective "first" and "second" sound speeds in explicit form. 

To see how this comes about, we proceed by evaluating the coefficients involved using 
the cool equation of state as specified by substitution of (|8^7| ) in ( |3.5|) . The lowest order 
coefficients can easily be worked out exactly as 

C = 3(c 2 p -v 2 )Q, C 00 = (3v 2 - cl)c 2 p Q , 
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C n = (3c 2 - f )Q , C 01 = -2vc 2 p Q , (A19) 

in terms of a factor 

Q = ^(cl-v 2 r i/3 c; 1/3 s s - 2/3 , (A20) 

which is unbounded, growing proportionally to G -2 , in the limit as — > 0. The coefficients 
of next order are bounded in this limit, but are functionally more complicated. Nevertheless, 
to the order of accuracy that is needed they are simply expressible in terms of the dilation 
coefficient $ 2 defined by ( |3.10|) as 



B = $ 2 , B 00 = $ 2 + 0{G 4 } , B u = $ 2 + O{0 4 } . (A21) 

c p 

To the required order of accuracy, the dilation coefficient itself (which effectively controles 
the dynamics of the zero temperature limit ||14|| ) will be given simply by 

71 

$ 2 = _ + O{0 4 } . (A22) 
A 4 

Those of the remaining coefficients that are not exactly zero are not just bounded but are 
characterised more strictly by 

A = 0{&} , A 00 = 0{&} , A n = 0{&} , 
A m = 0{&} , A 10 = 0{Q} , B 01 = 0{& 4 } . (A23) 

Thus they all tend to zero in the limit G — >• 0, and so will drop out as far as the present 
calculation is concerned. 

Since the combination Qs s 2//3 is bounded, it can be seen from (|A23| that the characteristic 



equation obtained as a condition of vanishing determinant for the combined system (|A14|) , 
( |A15|) , ( |A16| ), (|A17|) will be expressible to lowest order as a product of a pair of quadratic 



factors 

= b oo u 2 + c 2 {B u cos 2 9 + B sin 2 9) , ( A24) 

F = (C(C llU 2 + 2C 01 cos 6 u + C 00 cos 2 9) + (C 00 C n - C 01 2 ) sin 2 9 ) s s 4/3 
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in the form 

= 0{6 4 } , (A25) 

in which the left hand side is finite while the right hand side tends to zero as —>■ 0. At 
lowest order we thus obtain an effective decoupling into two factors one or other of which 
must vanish separately. The first alternative, T x = 0{O 4 }, will characterise ordinary "first" 
sound modes, and is expressible using the explicit expressions ( A21|) simply as 



u 2 ~c 2 p = 0{6 4 } . (A26) 

It can thus be seen that the "first" sound speed will be given, independently of the propa- 
gation angle 9, by 

u = ±c p + 0{Q 4 } , (A27) 

consistently with the original interpretation of c p as the ordinary sound speed in the zero 
temperature limit. The more interesting alternative, JF E = (9{B 4 } will characterise "second" 
sound modes, and is expressible using the formulae ( A19|) as 



(3c 2 - v 2 )u 2 - 4c 2 cos6vu-c p +c 2 p (l + 2 cos 2 6)v 2 = C{6 4 } . (A28) 
The "second" sound speed is thereby found to be given by 



2c 2 p vcos6±c pX /(c 2 p -v 2 )(3c 2 p - (1 + 2 cos 2 9)v 2 ) 
3c 2 — t> s 



u = - ry>1 ^lV — + Om . (A29) 



■ p 



It is to be remarked that when the two constituents are relatively at rest, i.e. when 
v = 0, this expression for the second sound speed reduces just to 

u = ±^ + O{0 4 } , (A30) 

whose form has been well known since the original work of Landau [§] in the Newtonian 
limit characterised by c 2 « c 2 . What is new here is the demonstration that this formula 
can be retained without any change for a relativistic superfluid in which the sound speed c p 
may be comparable with the light speed c. 
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